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Heavy-light hadrons are studied in a mass loaded flux tube model. The study indicates that 
the dynamics of mesons and baryons containing a c quark is well described by the mass loaded 
flux tube. The hypothesis of good diquark-antiquark degeneracy is found reasonable in heavy-light 
quark systems. The spectrum of charmed (D) and charmed strange (D a ) mesons is systematically 
computed. D and D s in ID multiplets are predicted to have lower masses in comparison with other 
theoretical predictions. The predicted masses of the l-(l 3 Di) and the3 _ (l 3 L> 3 ) D s agree well with 
those of recently observed D al {27QQ) ± and D s j(2860), respectively. 
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I. INTRODUCTION 

Hadrons spectrum can reveal the properties of the 
quark dynamics such as the color confinement. So far, 
a great progress has been made in the lattice QCD the- 
ory, but the quark dynamics in hadrons is not very clear 
and the hadrons spectrum can not be extracted from 
the QCD theory directly. The prediction of the hadrons 
masses has to be made in all kinds of models, and an ac- 
curate prediction would be a great challenge in hadrons 
spectroscopy. For heavy-light mesons, the spectrum has 
been systematically computed in the relativized quark 
model [lj], heavy quark symmetry theory Q, relativistic 
quark model [U, lattice QCD [5L chiral quark model @ 
and some other models p, 0, H, IE EE]- In these cal- 
culations, it is often difficult to predict the masses of 
higher orbital excited states. In many cases, the pre- 
dicted masses of the higher orbital excited states seems to 
be overestimated in comparison with experimental data. 

In hadrons containing more than two quarks or anti- 
quarks, two quarks or antiquarks may attract each other 
to make a diquark or anti-diquark cluster. The concept 
of diquark was put forth and was extensively studied in 
strong interactions [ll|, EE EE 

In terms of the diquark, a semi-classical mass loaded 
flux tube model [l9j] was recently exploited. In the model, 
a meson is considered a system with a massive quark m\ 
and a massive anti-quark m,2 connected by a flux tube 
(or relativistic string) with universal constant tension T 
rotating with angular momentum L. Similarly, a baryon 
is considered a system with a massive quark mi and a 
massive diquark to 2 connected by the flux tube. The 
flux tube is responsible for the color confinement. The 
mesons and the baryons are therefore described by the 
same dynamics in the same way. In addition, it is sup- 
posed that there is an approximate degeneracy between 
a good diquark in baryons and a relevant antiquark in 



mesons. 

Light mesons and baryons have been studied and clas- 
sified well in the model After the energy E and the 
angular momentum L of the system are written down 
with the dynamical parameters in the model, the E can 
be expressed in the L through some reductions. The gen- 
eral form of E is complicated, but the form is simple in 
some special cases. 

For the light quark systems, an approximate mass for- 
mula is given 



where T = £- 
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is the string tension, 
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m\ + m| with mi and mi are quark and anti- 
quark/diquark masses. 

The parameters for the light mesons and baryons in 
Eq. (TT|) were extracted from systematical analyzes of ex- 
isting data [IE H(| • The analyzes indicated that the pa- 
rameters for the light mesons match well with t he p aram- 
eters for the light baryons (see Table 7 in Ref. [19(). The 
dynamics (especially for large L) of light quark systems 
is well described by the mass loaded flux tube [Hj]. In 
the meantime, in order to account for an approximate de- 
generacy between the A baryons and the relevant mesons, 
the hypothesis of "good diquark-antiquark degeneracy" 
was proposed [l9j] ■ 

In heavy-light quark systems' case, an approximate 
mass formula was reduced also [IE ] 



E = M 
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where M is the heavy quark mass and m is the light 
quark/diquark mass, other parameters are indicated in 
Eq. (TTJ). The spin-orbit forces are ignored and L =/= in 
Eq. Q) and Eq. ©. 

However, the heavy-light quark systems have not 
been systematically analyzed except that some A c 
baryons (A c (2285), A c (2625), A c (2880)) were simply 
mentioned(with the relevant parameters M c = 1600 
MeV, m M = 180 MeV and a = 0.974 GeV 2 ) in Ref. [h|. 
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Many topics in the mass loaded flux tube model 
have not been studied in the heavy- light quark systems. 
Whether the dynamics of the mesons and the baryons can 
be described well by the flux tube has not been examined. 
The spectrum of charmed and charmed strange mesons 
has not been obtained. The hypothesis of "good diquark- 
antiquark degeneracy" , which holds in light quark sys- 
tems, has not been tested. 

In this article, the heavy-light quark systems are stud- 
ied in the mass loaded flux model with the inclusion of 
the spin-orbit interactions. The spectrum of the mesons 
containing one heavy c quark/ antiquark is systematically 
computed, and some possible interpretations of recently 
observed states are discussed. 



II. CHARMED AND CHARMED STRANGE 
MESONS 

In the conventional quark model, mesons may be 
marked by their quantum numbers n 2S+1 Lj, where n 
is the principle quantum number, S is the total spin, L is 
the orbital angular momentum, and J is the total angu- 
lar momentum. In most quark models, the interactions 
between the quark and the antiquark include the spin- 
independent confinement interaction, the spin-dependent 
interactions (spin-orbit interaction, color hyperfine inter- 
action) and some other interactions [l], Hi], [22|, [23| . The 
spin-orbit interaction consists of a color-magnetic piece 
and a Thomas-precession piece. The spin-orbit interac- 
tion is often considered the dominant one except for the 
confinement interaction, and is sometimes simplified an 
L ■ S coupling. This kind of spin-orbit interaction is em- 
ployed in our study, while other spin-dependent interac- 
tions such as the spin-spin interaction will be ignored. 

If the spin-orbit interaction was added to the energy 
E of the system from the beginning, the final relation 
between the E and the L would be much more compli- 
cated than Eq. @. As a good approximation, a term 
aL ■ S responsible for the spin-orbit interaction can be 
brought into Eq. @ phenomenologically. The parameter 
a is assumed a constant for the mesons having the same 
flavors(a depends mainly on the heavy flavor). It can be 
determined from the fit of experimental data. The study 
in this article indicates that the mass loaded flux tube 
with the inclusion of the aL ■ S coupling is very potential 
to produce the whole D and D s spectrum comparable to 
the experimental data. 

As well known, the heavy quark symmetry applies in 
the heavy- light mesons. In the heavy quark limit, the 
mass and spin sq of the heavy quark decouples. All the 
mesons properties arc determined by the light degrees 
of freedom. The spin-parity j p (total angular momen- 
tum j = Sq + I of light degrees of freedom) are good 
quantum numbers and are conserved in strong interac- 
tions. In heavy quark effective theory (HQET), the spin- 
dependent interactions depend on j. A natural way to 
account for the spin-orbit interaction in HQET is to in- 



clude the al ■ Sq coupling instead of the aL ■ S coupling. 
However, from our analysis of the experimental data, the 
D mesons spectrum is difficult to be reproduced with 
the simple inclusion of al ■ Sq in Eq. Besides, there 
may exist a spin-orbit inversion problem 24] in HQET. 
The heavy quark symmetry seems a little difficult to be 
accommodated in the present flux tube picture. This 
difficulty is left as an open question and is not studied 
here. 

The reason of the inclusion of the aL ■ S coupling can 
be realized in another way. The inclusion of the aL ■ S in 
Eq. ([2]) will result in a nought of hyperfine splitting (spin- 
triplet and spin-singlet splitting) of P-wave or D-wave 
multiplet, which is consistent with theories and experi- 
ments. These hyperfine splitting relations have already 
been predicted in many quark potential models. The 
hyperfine splitting relations hold very well in P-wave or 
D-wave multiplets of charmonium (even in IP multiplets 
of D mesons) [25j . 

Therefore, it is reasonable to extend the mass formula 
of the heavy-light quark systems to 

E = M+ \[^y + 2iKL _ im5 +a L ■ S (3) 

with 

f g_ J(J+l)-L(L + l)-S(S+l) 
L-S- - . 

With this formula in hand, we go ahead with the study 
of the heavy-light mesons. Firstly, we examine whether 
the hypothesis of "good diquark-antiquark degeneracy" 
is favored. For this purpose, the parameters M c = 1600 
MeV and a = 0.974 GeV 2 extracted from the charmed 
baryons [l^] are used as our inputs to compute the spec- 
trum of the D mesons. Under the hypothesis of good 
diquark-antiquark degeneracy, m UtC t = = 180 MeV. 

We obtained m(l 1 Pi) = 2.406 GeV for one IP D meson. 
This predicted mass agrees well with that of the exper- 
imentally observed -Di(2430)° (25|. Spectrum of other 
charmed mesons (L > 0) can be subsequently computed 
after a — 24.6 MeV has been fitted from other three IP D 
triplets. In terms of these parameters, m s — 320 MeV is 
determined from P/ S i(2536) ± and P/ S 2(2573) ± . The spec- 
trum of D s mesons can be systematically computed(the 
results are not given here for the reason mentioned be- 
low). 

The experimental spectrum of D and D s mesons can 
be well reproduced by the same group of parameters from 
the charmed baryons except that the predicted l 3 Po and 
l 1 Pi D s mesons are much heavier in comparison with 
possible experimental candidates. The fact that the spec- 
trum of the D, D s mesons and the charmed baryons is 
successfully obtained by the same formula and parame- 
ters indicates explicitly that the dynamics of the mesons 
and the baryons containing one heavy c quark is well de- 
scribed by the flux tube. It is found that the hypothesis 
of "good diquark-antiquark degeneracy" is favored in the 
heavy-light quark systems. 
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In Ref. the as are a little different for different 
kinds of light mesons and baryons, for which there are two 
reasons. One reason is that the <rs in the reference were 
extracted with spin-orbit interactions ignored. The other 
reason is that the string tension(the string is responsible 
for the dynamics) may be different for hadrons containing 
different flavors. Therefore, the as for mesons may be 
different from the as for baryons. In order to compute 
the spectrum of D and D s in a more reasonable way, 
the parameters of a and a have to be refitted from the 
confirmed D mesons. 

For a consistent study, the masses of the c quark 
and the light u, d quarks are regarded universal param- 
eters for the mesons and the baryons in our fitting pro- 
cesses. That is to say, the parameters m c = 1.6 GeV 
and m Ut d = 180 MeV extracted from the A c baryons are 
used as inputs to predict the spectrum of the charmed 
mesons. Other parameters a — 1.10 GeV 2 and a = 37.9 
MeV are extracted from the four IP charmed mesons 
candidates (to extract this two parameters, the minimum 
of mean square error of the mass of the four IP charmed 
mesons is applied). In terms of these parameters, it 
is straightforward to get the spectrum of IP and IP 
charmed mesons from Eq. (|3|). 

In experimental side, each observed state has a mass 
uncertainty. The mass uncertainties of observed states 
may result in some uncertainties to our predictions. If 
a mass uncertainty ±30 MeV is assumed for each IP 
charmed candidate, the a will have an uncertainty ±0.09 
GeV 2 . This assumed uncertainty may result in ±30 MeV, 
±44 MeV and ±54 MeV uncertainties to the masses of 
the IP, ID and IP charmed and charmed strange mul- 
tiplcts, respectively. 

Our results for the charmed mesons are obtained in Ta- 
ble 1. In the table, possible candidates for the D mesons 
of each state are displayed. For some states, there is 
no one to one correspondence between the j p and the 
n 2S+1 Lj notation. To compare our results with other 
theoretical predictions explicitly, we listed the results of 
two typical computations [J In Ref. [![, the nota- 
tion n 2S+1 Lj was used, and this notation is employed in 
our calculation. The calculation was performed in HQET 
and the notation j p was used in Ref. :5|. For simplicity, 
quantum numbers J p , j p and n 2S+1 Lj are all labeled 
in the table. A parenthesis is put for the j p when there 
is no one to one correspondence between the j p and the 
n 2S+1 Lj notation. A dash "-" is put in the entry where 
the corresponding state has not been computed in the two 
models. A "?" indicates that there is no observed candi- 
date corresponding to the assignment at present time. 

Our results for the IP states are comparable with 
those in Refs. [l], H| and experiments. For the ID states, 
our results are much lower in comparison with those in 
Refs. [H,i. This obvious difference may provide a way to 
examine whether the mass loaded flux tube model is rea- 
sonable. It may give people a hint to find an underlying 
dynamics of hadrons. 

In terms of the parameters a = 1.10 GeV 2 , m c = 1.6 
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TABLE I: Spectrum of D mesons(GeV) with parameters a = 
1.10 GeV 2 , m c = 1.6 GeV, m u , d = 180 MeV and a = 37.9 
MeV. 



GeV and a — 37.9 MeV extracted from the charmed 
mesons and baryons, the strange quark mass m s = 288 
MeV is determined from two IP charmed strange mesons: 
P s i(2536) ± and P s2 (2573) ± . The spectrum of the D s is 
subsequently computed and listed in Table 2. 

P* (2317) ± and P s i(2460) ± are two "exotic" states. 
They were firstly observed by BaBar [2^, [26[ and 
CLEO (2f| [2?J and were once interpreted as the + | 

and the 1 + j + D s mesons, respectively. However, there 
are different interpretations to them. One difficulty of the 
D s mesons interpretation is that they have lower masses 
in comparison with theoretically predicted masses. So 
far, this two states have not yet been pinned down def- 
initely. In our article, they are not used as inputs to 
determine the mass of the strange quark. The difficulty 
of the D s mesons interpretation is not yet solved in the 
mass loaded flux tube. P* (2317) ± and P s i(2460) ± are 
really difficult to be interpreted as the pure l 3 Po and 
1 Pi(l 3 Pi will mix with 1 1 P\) D s mesons. 

The situation of the D s mesons is similar to that of 
the D mesons. The predicted masses of the ID D s are 
much lower than those in Refs. p], H|. 

Recently, two new D s candidates were observed. 
P s j(2860) was first reported by BaBar [H in 

P sJ (2860) -> D°K+ , D+K° s 

with M = 2856.6 ± l.5(stat) ± 5.0{syst) and r = 48 ± 
7(stat)±10(syst) MeV. For its natural spin-parity: J p = 
0°, 1~, • • • , this state was explained as the first radial 
excitation of the P* (2317) or the 3~(1 3 P 3 ) [US, [13? 
]■ 

X(2690) was also reported by BaBar [28j], but the sig- 
nificance of the signal was not stated. 

P sJ (2700) was first observed by Belle [H in 

P+ -> D°D sJ -> D°D a K+ 
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TABLE II: Spectrum of D 3 mesons(GeV) with parameters 
a = 1.10 GeV 2 , m c = 1.6 GeV, m s = 288 MeV and a = 37.9 
MeV. 



with M = 2715 ± 11+^ and F = 115 ± 20 36 32 MeV. 
The mass and the decay width change a little in their 
published version [33[. For its J p = 1~, this state was 
interpreted as a mixture of the 2 3 S\ and the l 3 £>i [3(3] 
or the l-(l 3 L>i) HH. 

In these interpretations, one difficulty for the l 3 Di 
and the 1 3 D^ D s interpretations is that the masses of 
L> s i(2700) ± and D sJ (2860) are 100 200 MeV lower 
than the theoretical predictions. In the mass loaded flux 
tube model(Table 2), there is no difficulty to these in- 
terpretations at all. The predicted mass of the l 3 £>i 
D s is around 2714 ± 30 MeV and the predicted mass of 
the 1 3 L> 3 D s is around 2903 ± 44 MeV. As the masses 
and the decays modes considered, it is very possible that 
Z? s i(2700) ± and D sJ (2860) are the l 3 £>i and the 1 3 £> 3 
charmed strange mesons, respectively. 

III. CONCLUSIONS 

In summary, the mass loaded flux tube with the in- 
clusion of the spin-orbit interaction is studied. In heavy- 



light quark systems, the mesons and the baryons are well 
described by the mass loaded flux. Experimental data 
(spectrum) for the masses of the mesons and the baryons 
containing one heavy c quark are reproduced well by the 
same formula and the same parameters. It is found that 
the hypothesis of "good diquark-antiquark degeneracy" 
is a reasonable and consistent hypothesis in heavy-light 
quark systems. 

Our results indicate that 1^(2317)* and D s i(2460) ± 
are unlike the pure l 3 Pb and 1 Pi charmed strange 
mesons, respectively. 

Our predictions of masses of the II? D and D s are 
much lower in comparison with other theoretical pre- 
dictions. The predicted masses of the l~(l 3 Z?i) and 
the 3~(1 3 -D3) charmed strange mesons agree well with 
those of the recently observed D s i(2700) ± and D s ,/(2860) 
states, respectively. Other two ID charmed strange 
mesons around 2800 MeV are expected to exist, which 
is left for the confirmation of future experiments. 

Of course, many observed states are mixed states in the 
real world. Under mixing, how to interpret the observed 
states with the pure states is not clear, which deserves 
more study. 

The heavy-light quark systems containing one b quark 
have not been analyzed. The heavy quarkonium has not 
been explored either. The systems with radial excita- 
tion or excitation inside the string are not yet involved. 
How to extend the model to compute the spectrum of 
all kinds of mesons and baryons would be an interesting 
work. Furthermore, how to develop the model to describe 
the production and decay dynamics deserves further ex- 
ploration. It will be more important to find whether there 
is an underlying dynamics in the mass loaded flux tube 
model different from existing QCD inspired models. 

The mass loaded flux tube model is a semi-classical 
one, it will be interesting to study the mass loaded flux 
tube in a "fundamental" theory such as the string the- 
ory(some features such as the Regge trajectory behavior 
in the mass loaded flux tube model have already been 
obtained in the string theory). 
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